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Abstract. First detections of thermal water vapor ab- 
sorption lines have been made toward Orion IRc2 using 
the Short Wavelength Spectrometer (SWS) on board the 
Infrared Space Observatory (ISO). Grating spectra cover- 
ing wavelengths 25-45 /im yield 19 pure rotational lines, 
originating from energy levels 200-750 K above ground. 
Fabry-Perot spectra of 5 transitions resolve the line pro- 
files and reveal the H2O gas kinematics. The fact that all 
lines are seen in absorption is in striking contrast with 
data from the ISO Long Wavelength Spectrometer (LWS), 
where the H2O lines appear in emission. At least one line 
displays a P-Cygni type profile, which suggests that the 
water is located in an expanding shell centered on or near 
IRc2. The expansion velocity is 18 km s _1 , in agreement 
with the value inferred from H2O maser observations by 
Genzel et al. (1981). Because the continuum is intense 
and likely formed in or near the water-containing gas, the 
excitation of the observed transitions is dominated by ra- 
diative processes. A simple, generalised curve-of-growth 
method is presented and used to analyze the data. A mean 
excitation temperature of 72 K and a total H2O column 
density of 1.5 x 10 18 cm -2 are inferred, each with an esti- 
mated maximum uncertainty of 20%. Combined with the 
H2 column density derived from ISO observations of the 
pure rotational H2 lines, and an assumed temperature of 
200-350 K, the inferred H 2 abundance is 2-5xl0~ 4 in 
the warm shocked gas. This abundance is similar to that 
found recently by Harwit et al. (1998) toward Orion using 
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data from the LWS, but higher than that found for most 
other shocked regions by, for example, Liseau et al. (1996). 
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1. Introduction 

Water is one of the prime species for probing the inter- 
action between young stars and their surroundings, both 
in terms of its abundance and its peculiar excitation. In 
high temperature gas, appropriate to shocks and hot cores 
for instance, all of the oxygen not locked up in CO is pre- 
dicted to be driven into H2O at temperatures above 230 
K by the O + H 2 -> OH + H and OH + H 2 -> H 2 + H 
reactions, resulting in very bright H 2 lines (e.g., Hollen- 
bach & McKee 1979, Neufeld & Melnick 1987, Kaufman & 
Neufeld 1996, Charnley 1997). In contrast, the H 2 abun- 
dance may be at least two orders of magnitude lower in 
surrounding colder gas (e.g., Zmuidzinas et al. 1995). In 
addition to collisions, the H 2 excitation and line profiles 
can be strongly affected by mid- and far-infrared radiation 
from warm dust (e.g., Phillips et al. 1978, Takahashi et al. 
1983, 1985), providing detailed information on the physi- 
cal parameters of the gas and its location with respect to 
the radiation sources. 

Interstellar water has been difficult to detect, apart 
from its presence as an ice mantle on dust grains or 
through its maser emission, owing to the severe telluric ab- 
sorption encountered at Earth-based observatories. Never- 
theless, over the last twenty years, many searches for lines 
of gas-phase H 2 and its isotopomers have been made 
from the ground and airborne altitudes, in particular to- 
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ward Orion (e.g., Waters et al. 1980, Phillips et al. 1978, 
Jacq et al. 1990, Wannier et al. 1991, Cernicharo et al. 
1994, Zmuidzinas et al. 1995, Tauber et al. 1996, Timmcr- 
mann et al. 1996, Gcnshcimcr et al. 1996, Cernicharo et al. 
1999a, b). However, due to the choice of line, wavelength 
and beam size, different components are often probed, 
whilst many ground-based observations refer to masing 
lines, for which sophisticated shock and maser models are 
required in order to extract physical parameters. 

One of the major aims of the Infrared Space Observa- 
tory (ISO) mission has been the routine measurement of 
thermal gas-phase water lines, and their use as diagnostics 
of the chemical and physical conditions within molecular 
clouds. Far-infrared pure rotational H 2 emission lines 
in the 50-200 /im wavelength range have been detected 
with the Long Wavelength Spectrometer (LWS) of ISO in 
a number of star-forming regions (e.g., Liseau et al. 1996, 
Ceccarelli et al. 1998), including Sgr B2 (Cernicharo et al. 
1997a) and Orion (Cernicharo et al. 1997b, 1999a; Har- 
wit et al. 1998). Additionally, van Dishoeck & Hclmich 
(1996), van Dishoeck (1998) and Dartois et al. (1998) 
have observed absorption around 6 /im in the z^ 2 = l 
band toward a number of deeply embedded, massive young 
stars. Typical H 2 abundances of 10~ 5 with respect to H 2 
have been derived from these data. Similar observations 
have recently been reported for Orion BN/IRc2 by van 
Dishoeck et al. (1998) and Gonzalez- Alfonso et al. (1998), 
although in this case emission is also detected. Wright et 
al. (1997) have however shown that the detection of the 
corresponding pure rotational lines at ~ 30-200 /<m in 
most sources observed at 6 /im is still difficult. The ob- 
servations of Orion-IRc2 presented here form a notable 
exception. 

Many of the earlier searches for H2O lines have 
been performed toward the BN/IRc2 complex of infrared 
sources in Orion, because of the extraordinary brightness 
of many atomic and molecular lines in this region com- 
pared with other clouds (e.g., Genzel & Stutzki 1989, 
Blake 1997). See van Dishoeck et al. (1998) and references 
therein for a detailed description of the geometry and the 
diverse range of physical phenomena in this region. In this 
paper we present the first detection of numerous pure ro- 
tational water lines in absorption toward IRc2 in the 25- 
45 /im interval with the Short Wavelength Spectrometer 
(SWS) (de Graauw et al. 1996) on board ISO. Some of 
these lines have been velocity resolved with the Fabry- 
Pcrot, enabling direct information on the location of the 
absorbing gas to be inferred. These data complement the 
earlier ground-based data, as well as observations of pure 
rotational lines with the LWS obtained by Cernicharo et 
al. (1997b, 1999a) and Harwit et al. (1998) in a much 
larger beam. 

Our interpretation of the ISO spectra of Orion IRc2 
suggests that the 25-45 /tm H 2 spectrum originates in 
a region where the intrinsically strong lines couple effi- 
ciently to an intense continuum. In principle, the forma- 



tion of such a spectrum should be described for a stratified 
atmosphere in which lines and continuum are treated self- 
consistently. We show here that the observed features of 
the spectrum can be described well by a simple "gener- 
alised curve-of-growth" , which includes the effects of cou- 
pling to a strong continuum. In the following we describe 
our observations and present our results, followed by a 
discussion of the location of the absorbing water vapour, 
its excitation and finally its abundance. 

2. Observations and data reduction 

2.1. Observations 

A complete grating scan from 2.4-45.3 /im using the 
maximum spectral resolution SWS06 observing mode was 
made on September 6 1997 (revolution 660) centered at 
a(2000) = 05 /l 35 m 14.2 s , <5(2000) = -05°22'31.5", which 
is about 3" W and 1" S of the IRc2 position quoted by 
Downes et al. (1981). The long axis of the ISO-SWS aper- 
ture was oriented at 172.81° E of N. The aperture size 
varies from 14" x 20" at 2.4-12 /im, to 14" x 27" at 12- 
27.5 /jm, 20" x 27" at 27.5-29 /im and 20" x 33" at 29-45.2 
/im. A further SWS06 scan was made from 26.3 to 45.2 /im 
in revolution 826 (February 19 1998) with the long axis 
oriented at 165.79° E of N. The SWS beam includes both 
IRc2 and BN, but not "peak 1" or "peak 2" of shocked 
H 2 (Beckwith et al. 1978). 

The resolving power of the ISO-SWS grating, A/AA, 
varies from about 1000-2500, implying that the observed 
line profiles are not resolved and that little velocity in- 
formation can be obtained to help disentangle the various 
components which are known to exist toward IRc2. For 
this reason, follow-up observations of a selection of H 2 
lines were obtained using the Fabry-Perot SWS07 observ- 
ing mode in revolutions 823 and 831 (February 16 and 24 
1998). Five 250 km s^ 1 scans were made across each line. 
The resolving power is of order 30 000, corresponding to 
a velocity resolution of 10 km s _1 , and the aperture long 
axis was oriented at 164.14° and 168.46° E of N on the 
two dates. The aperture size is 10" x 39" and 17" x 40" for 
wavelengths below and above 26 /im respectively. During 
the revolution 831 observation the wavelength interval be- 
tween 5.3 and 7.0 /im was simultaneously scanned by the 
SWS grating. Within this region are several H 2 lines (0-0 
S(7), S(6) and S(5)), and H 2 solid state ice and H 2 gas 
phase i/ 2 =l-0 features. 

2.2. Data reduction 

Data reduction was carried out on the Standard Processed 
Data file from the Off Line Processing system (OLPv6.1.1 
for revolution 660 and OLPv6.3.2 for revolutions 823, 826 
and 831), using standard routines within the SWS Inter- 
active Analysis package up to the Auto Analysis Result 
(AAR) stage. The dark current subtraction was performed 
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interactively for the grating data, but not for the FP data. 
The most up-to-date calibration files available for wave- 
length, flux and relative spectral responsivity were used. 
Even so, there is at least one instrumental artefact in our 
grating spectra, at 33 /im, due to structure in the relative 
spectral response calibration file (RSRF). 

The high continuum flux encountered in band 4 of the 
SWS grating (29-45 /im), and the finite response time 
of the detectors, led to significant memory and response 
effects near the band edges, at positions of important wa- 
ter and OH lines, and so necessitated a non-standard up- 
down scan correction. For the SWS an "up" ("down") 
scan proceeds toward increasing (decreasing) grating po- 
sitions but decreasing (increasing) wavelengths. Since the 
up scan precedes the down scan, in high flux cases it is 
usually the down scan which is deemed to have the more 
reliable spectral shape, since by that time the instrument 
has settled down. Indeed, this was the case in the long- 
wavelcngth portion (41-45 /im) of our spectra, so that the 
up scan was corrected to approximately the same shape as 
the down scan. On the other hand, the short-wavelength 
portion (29-33 /im) of the down scan of our spectra had 
a peculiar concave shape, which did not match cleanly (in 
terms of shape) to band 3e, which in turn matched cleanly 
with band 3d. This was true for both the revolution 660 
and 826 data sets. In this case therefore the down scan 
was corrected to have a similar shape to the up scan. The 
up-down correction obviously entailed splitting the band 
4 data into two segments (29-37 /im and 37-45 /im), but 
this produced a more reliable spectral shape at both ends, 
e.g. allowing the continuum on either side of the 45.1 /im 
water line to be determined, and the 28.940 and 28.914 
/mi OH and H 2 44 — 3 13 lines to be partially resolved. 

Following the AAR stage, further data processing, such 
as fiatfielding, sigma clipping and rcbinning, was per- 
formed using software in the SWS IA package. For in- 
stance, for both the grating and FP cases the individual 
detector scans have been brought to a common flux level 
using an n th degree polynomial (n = 0, 1, 2 or 3) fit to a 
reference spectrum which may either be the median of the 
down scans (e.g. grating band 3), the median of both the 
up and down scans (e.g. the up-down corrected grating 
band 4), or the mean of all scans (e.g. the Fabry-Perot 
data). The data have subsequently been sigma clipped 
such that any points lying more than 3 sigma outside of 
the average of all data within a bin of width equal to 
the theoretical spectral resolution have been discarded. 
The final step in producing a spectrum involves rebinning 
the much over-sampled data using the mean within a bin 
which is again equal to the theoretical resolution. 

In the case of the 25.94 /im FP scans it was observed 
that glitches and their associated tails, which had not been 
flagged as bad data by the pipeline, contributed significant 
noise to the final spectrum. In this instance an interactive 
glitch and tail correction was performed before fiatfielding, 
sigma clipping and rebinning. As noted by Schaeidt et al. 



(1996) and Heras (1997) the photometric accuracies of the 
SWS grating band 4 (where most of our water detections 
occur) and Fabry-Perot are 30% and 40% respectively. 
However, our observed line equivalent widths, being ratios 
of the line area over adjacent continuum, will be more 
accurate than this. The accuracy of the grating wavelength 
calibration is 1/10-1/5 of a resolution clement, whilst that 
of the Fabry-Perot is of order 10~ 4 /im, or about 1 km s _1 
(Valentijn et al. 1996; Feuchtgruber et al. 1997). 

3. Results 

In Fig. 1 the SWS grating spectrum in the 25-45 /tm range 
is presented for the revolution 660 data, with the position 
and identification of features marked. The data consist 
of spectral segments from bands 3d, 3e and 4, for which 
"jumps" in the flux were observed due to the different 
aperture sizes. However, for display purposes small shifts, 
of order a few thousand Jy (i.e. a few percent), have been 
applied to produce a continuous spectrum. Clearly marked 
in Fig. 1 are 19 absorption lines arising from pure rota- 
tional levels in both the ortho- (o) and para- (p) species 
of H 2 0. In the few cases where the plot scale makes the 
line difficult to see in Fig. 1, we present close-ups in Fig. 2. 
The lines in this wavelength region originate from highly- 
excited states, with excitation energies of the lower lev- 
els hcEz/k = 200 — 750 K, and the upper levels up to 
1200 K. The strongest absorptions come from levels be- 
longing to the backbone of the energy ladder. Among the 
highest levels seen is the upper state of the 7 2 5 — 616 tran- 
sition at 29.8367 /tm, involving the upper level (616) of the 
well-known 22 GHz H 2 maser line (Cheung et al. 1969). 
Indeed, many of the levels from which absorption is seen 
are also upper levels of either known or suspected maser 
transitions (Ncufcld & Melnick 1991). It is shown below 
based on velocity considerations that our lines probably 
arise in the same region where the masers are situated, 
indicating that our observations may form a useful data 
set for studying the maser pump mechanism. 

A sub-set of 5 of the 19 lines was chosen to be ob- 
served with the Fabry-Perot. These spectra are shown in 
Fig. 3. In one case (4 32 — 3o3 40.69 /im) and possibly a 
second (541 — 4 32 43.89 /im), two velocity components are 
detected, one in absorption and the other in emission. A 
modulation in the baselines, due to the tracking of the 
Fabry-Perot transmission function across the region of the 
peak of the grating instrumental profile, has been mod- 
elled in terms of two parameters, the grating spectral res- 
olution and the shift between the two response functions. 
The spectra shown in Fig. 3 have been corrected for this 
tracking pattern. Also shown in Fig. 3 is the grating spec- 
trum from 5.3-7.0 /tm obtained from these new data. It is 
consistent with that already observed by van Dishoeck et 
al. (1998) and Gonzalez-Alfonso et al. (1998). 

Table 1 presents an overview of our observations for the 
Revolution 660 data, including the rest wavelength, upper 



4 



CM. Wright et al.: Water absorption toward Orion IRc2 



5.0 

4.5 

~~> 4.0 
"o 

3.5 

X 

^ 3.0 
2.5 
6.5 

~ 6.0 
—> 

*o 

5.5 
_^ 

Li_ 

5.0 
8.0 

o 

° 7.0 

_3 

Li_ 

6.5 

9.5 

? 9.0 
o 

x 8.5 

ni 

8.0 







! ^2® 63 3 — 5 24 



40 



41 



42 43 44 

Wavelength (/Ltm) 



45 



46 



Fig. 1. ISO-SWS 25-45 /im grating spectrum toward Orion IRc2. The principal absorption and emission features are 
indicated. The data within bands 3d, 3e and 4 have been shifted so that they join continuously at the band edges. The 
feature marked RSRF is an instrumental artefact due to structure in the relative spectral response calibration file. 



and lower levels, and equivalent widths W = cW\/\ in 
units of Doppler velocity (km s _1 ), where W\ is the line 
equivalent width in wavelength units. The superscripts on 
W refer to the observed and calculated equivalent widths, 
the latter derived from a generalised curve-of-growth anal- 
ysis to be described in Sect. 4.2. Table 1 also contains 
the results for three pure rotational OH lines detected 
in absorption in our spectra. All derived quantities were 
calculated from unshifted spectral segments. Figure 4(a) 
presents our results graphically, in the form of a standard 



population diagram (Goldsmith & Langer 1999), where 
the column density in the lower level has been determined 
using the formula N t = (8ncW\ge)/(\ 4 A u eg u ), where g e 
and g u are the statistical weights of the lower and upper 
levels, A U £ is the Einstein A-cocfficient in s _1 and W\ is 
the equivalent width. This formula is only strictly valid 
for optically thin lines, and under the assumption that 
the upper levels are not significantly populated. The er- 
ror bars represent a 30% uncertainty, so that the scatter in 
the data is clearly larger than any systematic or statistical 
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Fig. 2. Close-ups of several of the H2O absorption lines shown in Figure 1. The rest wavelength of each line is indicated 
by a tick mark. The P-H2O 642-533 39.3992 /im line is visible as a shoulder on the deeper 0-H2O 550-441 39.3749 /im 
line. 



uncertainties. Although there is a clear trend in the data, 
this scatter already hints that a more detailed analysis is 
required to obtain a reliable total H 2 column density 
Such an analysis, in the form of a generalised curve-of- 
growth, is presented in Sect. 4.2. 

Almost all lines listed in Table 1 at 27.5-45.2 /im are 
visible in both the Revolution 660 and 826 spectra, and 
their measured equivalent widths agree within 30% in the 
majority of cases. See Fig. 4(a), where both data sets are 
presented graphically. We present only the Revolution 660 
data in Table 1 because the data from that revolution are 
of a better quality, and also only 15 of the 19 lines were 
observed in Revolution 826 due to the shorter wavelength 
scan range. The fact that the pure rotational H2O and 
OH lines observed with the SWS occur in absorption is in 
strong contrast with the lines seen with the LWS, which 
are mostly in emission (e.g., Harwit et al. 1998). Since 
the LWS lines arise in levels with similar excitation ener- 
gies at wavelengths where the continuum flux is compara- 
ble to that in the SWS range, this suggests that most of 



the mid- infrared continuum and H2O absorption originate 
in a small region of size comparable to the SWS beam. 
In the almost order-of-magnitude larger area of the LWS 
beam, even more of the absorption is probably filled in by 
large-scale emission from more extended, lower excitation 
shocked gas. The qualitative change in the character of 
the water spectrum is thus probably a result of effects of 
both geometry (beam-filling) and radiative transfer. 

The small differences in equivalent widths found from 
the grating and Fabry-Perot spectra are due to uncertain- 
ties in the calibration and varying beam shapes, which 
may couple differently to the complex source structure. 
The continuum radiation at 25-45 /im is produced by ex- 
tended dust emission over a ~ 30" region centered near 
IRc2, in contrast with the situation at ~ 6 /im where BN 
dominates (Wynn- Williams et al. 1984). Since H2O emis- 
sion is observed on an arcmin scale with the LWS (e.g. 
Cernicharo et al. 1999a), it is likely that the H 2 covers 
all of the continuum, at least for the lower-excitation lines. 
However, the absorption of the lower-lying lines may be 
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Fig. 3. ISO-SWS Fabry-Perot spectra toward IRc2 of a selection of H2O absorption lines. The final panel displays 
the 5.3-7.0 /*m grating spectrum of the H2O f2 = l-0 band taken simultaneously with the revolution 831 Fabry-Perot 
spectra. 



partially filled in by emission within the SWS beam. As 
discussed below, these effects must be taken into account 
explicitly in an analysis of the equivalent widths. 

4. Analysis and discussion 

4-1. Spatial origin of the water 

It is important to establish where the observed water is 
located along the line-of-sight towards IRc2. Possible can- 
didates are the quiescent ridge, the shocked (low and/or 
high velocity) plateau gas, the hot core, or a combina- 
tion of these components. Fortunately, our Fabry-Perot 
observations can aid in resolving the ambiguity. For the 3 
pure absorption cases, the wlsr is ~ — 8±3 km s _1 . This 
blucshift is similar to velocities found for molecules such as 
CO, C2H2 and OCS of —3 to —18 km s _1 in absorption by 
Scoville ct al. (1983) and Evans ct al. (1991), which they 
interpret as arising from the low-velocity plateau gas, rep- 
resenting a shell expanding about a point close to IRc2. 



The most likely candidate for the location of our H2O 
gas is therefore the shocked, low-velocity plateau gas. Our 
lines are also resolved with observed full-width at half- 
peak AV^ ~ 30 km s . Assuming that the intrinsic source 
profile is Gaussian and that the FP instrumental profile 
is an ideal Airy profile with FWHM of 10 km s _1 , the 
intrinsic line widths are of order 23 km s . This broad 
lincwidth excludes most other components, such as the 
quiescent ridge and hot core, which have AV < 20 km 
s _1 (Blake et al. 1987). H2O differs in this respect from 
the other molecules seen in absorption, which have AV of 
order only 3-10 km s _1 (Scoville et al. 1983, Evans et al. 
1991). 

The peak of the emission component of the 40.69 ^m 
line occurs at «lsr=+7 to ~ +12 km s _1 . This velocity 
range is similar to that found by Cernicharo et al. (1994) 
for the 31,3 — 2 2 ,n 183 GHz line, as well as for a selec- 
tion of H3, 8 lines by Zmuidzinas et al. (1995), and is in 
fact close to vlsr ~ 9 km s~ x appropriate for the dense, 
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a Orion shock model of Kaufman & Neufeld (1996), using lower oxygen and H2O abundances of 3.16 x 10~ 4 and 3.5 x 10~ 4 
respectively, cf. Harwit et al. (1998), and assuming a shock velocity of 37 km s _1 . 

b The net equivalent width, including both emission and absorption components. Individually, the observed absorption com- 
ponent is 11.07 km s _1 , whilst the emission component flux is 7.8xl0 -18 W cm -2 . 

c The net equivalent width, including both emission and absorption components. Individually, the observed absorption com- 
ponent is 6.11 km s _1 , whilst the emission component flux is 2.3xl0~ 18 W cm~ 2 . 



quiescent ridge in which IRc2 is embedded. The emission 
component is also resolved, with an observed linewidth of 
AV w 40 km s _1 , consistent with the intrinsic widths of 
~ 35 km s _1 of the water lines detected in emission at 
longer wavelengths with the ISO-LWS by Harwit et al. 
(1998). 

In the comparisons presented above the absorbing and 
emitting components have been considered separately. 
However, the profile as a whole for the 40.69 /xm line, 
and to some extent also the 43.89 /zm line, is very remi- 
niscent of a P Cygni-type profile, and is similar to those 
observed for several far-infrared OH lines toward Orion by 
Betz & Boreiko (1989) and Melnick et al. (1990). In fact, 
the total blue-shift of the absorption component for all 



5 lines, i.e. with respect to the rest velocity of the cloud 
vlsr ~ 9 km s _1 , is very close to the expansion velocity of 
the "low- velocity flow" , (or plateau, or "expanding dough- 
nut") of 18±2 km s" 1 quoted by Genzel et al. (1981) from 
proper motion studies of the 22 GHz H 2 maser emis- 
sion. Together with the fact that the emission component 
of the 40.69 /xm line is near the rest velocity of the cloud, 
this is precisely the definition of a P-Cygni profile and 
establishes almost without a doubt that the observed wa- 
ter arises from an outflow centered near IRc2, most likely 
from source "I" (Menten & Reid 1995). As noted by Gen- 
zel et al. (1981) this outflow, or low velocity plateau, can 
be traced from ~ 2 x 10 17 cm (30") to within 10 14 cm 
from its dynamical center. 
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Fig. 4. (a) Level populations of H2O derived from the ab- 
sorption lines and using the optically thin formula Ne = 
(8ncW\g£)/(X 4: A u ig u ). The solid circles represent the rev- 
olution 660 grating data, the asterisks the revolution 826 
grating data, and the open diamonds the revolution 823 
and 831 Fabry-Perot data. The dotted line represents a 
best fit 1 st order polynomial to the revolution 660 data, 
and implies an excitation temperature of 71.6 K. A sim- 
ilar fit to the revolution 826 data yields 75.9 K, whilst 
the Fabry-Perot absorption (emission) data yield 62.9 K 
(64.4 K). (b) Level populations of H 2 derived from the 
generalised curve-of-growth (GCOG) analysis of the revo- 
lution 660 data. The best fit GCOG analysis for the pop- 
ulations results in a distribution that is consistent with a 
thermal distribution at a constant excitation temperature 
of 72.5 K. The dotted line indicates a best fit 1 st order 
polynomial to the points with Eg < 600 K. 



We note that Takahashi et al. (1985) predicted that 
many H 2 lines from a linearly expanding spherical cloud 
would exhibit P-Cygni type profiles, although the only 
case in common with our observations, the 40.69 /im line, 
is not P-Cygni shaped in any of their models. However, 
they do not include shock emission in their calculations. 
Further, there is a difference between a true P-Cygni pro- 
file and the P-Cygni type profiles observed here and pre- 
dicted by Takahashi et al. A true P-Cygni profile results 
from an expanding stellar atmosphere with a central lu- 



minosity source, whilst the P-Cygni type profiles of far- 
infrared water lines are caused by the presence of warm 
dust throughout the expanding cloud. Further details on 
water line profiles in regions of intense dust emission may 
be found in Takahashi et al. (1985) or Doty & Neufcld 
(1997). 

4-2. Excitation of the absorbing water vapour 

4.2.1. Generalised curve-of-growth method 

The ISO spectra of Orion highlight the need for a non- 
traditional analysis of interstellar spectra. At shorter 
wavelengths (in the visible and ultraviolet) , molecular ab- 
sorption lines arise in gas that is physically very cold in 
relation to the radiation temperature of the background 
star, so that the absorption can be described with a classi- 
cal curve-of-growth analysis and the effects of stimulated 
emission can be neglected. At longer wavelengths (radio), 
an isolated molecular region typically exhibits collisionally 
excited emission lines on top of a weak continuum so that 
the most significant radiative coupling is to the cosmic 
background radiation at T c ^ r = 2.728 K. In the present 
case, we observe a molecular spectrum that is formed in 
the presence of a very strong continuum, i.e., where the 
radiation temperature of the continuum may be compara- 
ble to the excitation energies of the states involved. This 
implies both that the classical curve-of-growth analysis of 
cold absorbers detached from the background continuum 
is invalid and that the standard analysis of collisional exci- 
tation in competition with a very cold, dilute background 
is inadequate. Indeed, it is the near equality of continuum 
brightness temperature T ra< x(v) and molecular excitation 
temperature T ox that explains at least partly why the wa- 
ter spectrum of Orion IRc2 is predominantly in absorp- 
tion at A Si 50 /im and in emission at longer wavelengths. 
In detail, such spectra should be modelled as an extended 
atmosphere in which lines and continuum are treated con- 
sistently. We will outline briefly how it is possible to ex- 
tract important information on excitation and abundance 
through a simpler analysis of the unresolved lines based 
on a generalised curve-of-growth. Details of this method 
will be discussed elsewhere (Black, in preparation). 

Consider the situation in which a fractional area a c 
of a continuum source of intensity I c {v) is obscured by a 
column of gas in the foreground. The observed spectrum 
in the vicinity of a transition u <-> I is given by the sum 
of the unattenuated continuum (first term), the absorbed 
continuum (second term), and the emission of the cloud 
itself (third term): 

^obs(^) = (1 - a c )I c (v) 
+ a c I c {v)e- T ^ 

+ B v {T ut )(l-e- T ^) , (1) 

where t u i is the optical depth and B v (T u g) is the Planck 
function evaluated for an excitation temperature T u £. 
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This temperature is defined by the column densities 
of molecules in the upper (N u ) and lower (Ng) states 
through: 

Nu _ 9u / 
N e ~ 



■ exp 



( kT„i ) 



(2) 



ge v kT u e , 
The optical depth function for a single line can be written 
(e.g. Rybicki & Lightman 1979, Wannier et al. 1991) 

_ 7 A u£ Ne q u (_ , hv 



Tut (y) = 3.738 x 10" 



'ul 



AV ge 



H 

if v 



1 — exp( 



kT. 



))#") 



(3) 

where Ne is in cm 2 , i> u e — v u i/c is the wavcnumber of 
the line in cm -1 , and the line shape function is a gaussian 
of full- width at half-maximum AV in km s _1 , normalized 
so that <t>(v u i) = 1-0 • The term in parentheses is the 
correction for stimulated emission. 

The generalised equivalent width, in frequency units, 
is expressed as 



w: 



obs 



/ 



i c {y) - lohsW) 



dv 



w; 



obs 



J [(a. 



B v (T u e) 

W) 

■ 1 and T, 



dv 



(4a) 



(46) 



In the limit as a c ~ 1 and T u e << hv/k, the equivalent 
width approaches the classical form of cold absorbers de- 
tached from the background continuum. By convention, 
the equivalent width is positive for net absorption and 
negative for emission. Emission can arise either when T u e 
is large enough that B v (T u p) / I c {v) > a c or when there 
is a population inversion (N u /Ne > g u /gi) so that r u e is 
negative, as in a maser. 

For a single unresolved line of equivalent width W° hs , 
there is not a unique solution, because 4 parameters must 
be constrained: a c , N u , Ne, and AV. In the classic prob- 
lem, in which a continuum point source of intensity is 
located behind, and completely obscured by, a column of 
cold foreground gas, and there is negligible population in 
the upper level of an absorption transition, then a c = 1 
and N u /N e ~ respectively. In that case, two lines arising 
in a common level suffice to determine Ne and AV. 

4.2.2. Results of the generalised curve-of-growth analysis 

We have analyzed the equivalent widths using Eq. (4b) . If 
we take the continuum intensity to be approximately con- 
stant over the line profile, then the term a c —B v (T u e) / I c {v) 
can be taken outside the integral in Eq. (4b). We adopt an 
empirical continuum intensity I c {v) = /° bs /^beam, where 
/° bs is the flux density observed in the ISO-SWS aperture 
of solid angle fibcam- When I c (v) is fixed in this way, the 
equivalent width of a line depends on the parameters JV>, 
T u e, AV, and a c . We adopt an intrinsic width AV = 23 
km s -1 and assume that a c has the same value for all lines 
in the following analysis. Thus each line is characterized 
by two parameters Ne and T u e (or Ne. and N u ). The solu- 
tion is still underdetermined because there are 29 distinct 



energy states involved in the 19 observed transitions. For- 
tunately, there are several cases in which some of the same 
states are involved in two or more transitions. As a first 
approximation, it is assumed that a single value of the ex- 
citation temperature T u e = T eXi o applies to all transitions 
with lower-state energies Ee < 275 cm -1 . The computed 
equivalent widths (Eq. 4b) agree with the observed values 
within 33% for all but the weakest one of these lines, when 
a c = 1, T OXj0 = 72.5 K, and N total = 1.5 x 10 18 cm" 2 . This 
solution is not very sensitive to the value of the continuum 
covering factor: for a c = 0.9, we obtain T ox ,o = 71.0 K and 
Ntotai = 1-6 x 10 18 cm~ 2 , while for a c = 0.8, the corre- 
sponding results are 69.1 K and 1.7 x 10 18 cm~ 2 . 

With the populations of the lowest energy states (Ee < 
275 cm -1 ) fixed at T e ^. = 72.5 K, the column densities 
in higher states are adjusted until the standard deviation 
of the mean W ohs — W GCOG is minimized for the entire 
set of lines. In this solution, the observed W ohs and cal- 
culated W GCOG equivalent widths all agree within 16%, 
with the single exception of the 441 — 4 i4 line at 37.98 
fim. The calculated equivalent widths are listed in Table 
1. Note that the adjustment also improves the fit for the 
lines arising in the low-lying states because W GGOG is 
sensitive to column densities in both the upper and lower 
states. The result is a population distribution that is indis- 
tinguishable from a thermal distribution at 72.5 K for all 
states with hcEi/k ^ 700 K. A few of the observed levels 
at higher energies (hcEi/k > 700 K) are well constrained 
because they are the upper states of more than one tran- 
sition, and their populations start to clearly show effects 
of subthermal excitation, T u e « 55 K. Even though the 
adopted linewidth is large, the derived line-center optical 
depths, listed in Table 1, reach values as high as r = 12.8 
and 9.7 in the 45.11 and 40.69 /tm lines, respectively. The 
smallest values of optical depth are To « 0.06. Finally, or- 
tho and para states of H 2 are consistent with the same 
excitation temperature and with an ortho/para ratio of 3. 

The results of our generalised curve-of-growth analysis 
are shown graphically in Fig. 4(b), in the form of a stan- 
dard population diagram. The much reduced scatter in the 
data is clearly a marked improvement over the optically 
thin approach assuming no stimulated emission used in ob- 
taining Fig. 4(a). It is however interesting that the inferred 
excitation temperatures from the two approaches are sim- 
ilar. The excitation temperature of ~ 72 K is similar to 
the colour temperature of the dust between 20 and 100 /im 
inferred by Werner et al. (1976), which confirms the ex- 
pectation that pumping by infrared continuum radiation 
from dust plays an important role in the H2O excitation. 
In this case, there is no direct constraint on the kinetic 
temperature of the H 2 0-containing gas. Subthermal exci- 
tation by collisions at densities below the critical densities 
of the observed levels (~ 10 9 cm~ 3 ) can also result in ex- 
citation temperatures of order 50-100 K, but this merely 
re-affirms that radiative processes probably dominate the 
excitation. We also note that the excitation temperature 
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inferred from the emission components of the 40.69 and 
43.89 ^m lines is <~ 64 K, also implying a large contribu- 
tion by radiative processes to the populations. Moreover, 
the analysis suggests that T u n of the absorption lines is so 
close to the continuum radiation temperature in this wave- 
length range, that it is easy to see why the H 2 spectrum 
should go into emission at A > 50 /xm as the opacity and 
radiation temperature of the continuum decrease with in- 
creasing wavelength. 

Because the lowest state involved in the 25-45 /zm 
spectrum lies at hcEi/k = 197 K, we cannot determine 
the column densities in the most populated states di- 
rectly from these observations. However, if we assume 
that the observed excitation pattern applies to the unob- 
served lower states, then we infer a total column density 
of iV(H 2 0) = 1.5 x 10 18 cm" 2 averaged over the ISO-SWS 
beam toward Orion IRc2. 

The net absorption by water in Orion is sensitive to 
the competition between the molecular excitation and the 
continuum brightness. We have assumed that the true con- 
tinuum intensity in Eq. (4) is equal to the mean surface 
brightness observed in the spectrum (flux averaged over 
the aperture). In general this need not be true, and the 
analysis can include an additional correction factor for the 
coupling to the continuum. Tests show that if the contin- 
uum intensity is thus changed by a factor of two in either 
direction, the derived total column density and mean ex- 
citation temperature change by less than 20%, which we 
take as our maximum uncertainty. 

4.2.3. Shock models 

We have also compared our observed line equivalent 
widths with those computed for the Kaufman & Neufcld 
(1996) C-shock in Orion with a pre-shock density of 10 5 
cm -3 and shock velocity of 37 km s -1 (Kaufman & 
Neufcld, private communication) (see Table 1). The wa- 
ter abundance (with respect to H 2 ) used in this shock 
calculation was 3.5xl0~ 4 , based on the O and C abun- 
dances of Cardelli et al. (1996), with all the carbon locked 
up in CO and H2O accounting for all the remaining oxy- 
gen. Although this model refers to "peak 1" rather than 
the IRc2 region, it reproduces the equivalent widths quite 
well within a factor of a few, except for some of the higher- 
lying lines (e.g., 27.03, 39.37, 40.76, 44.19 /an). Most of 
these lines are likely to be strongly affected by radiative 
excitation, which was not included in their model. 

We note that the shock velocity used is 37 km s _1 , dif- 
ferent from the expansion velocity inferred from our FP 
observations of 18 km s _1 . Indeed, our modelling of our 
H 2 emission, using the Kaufman & Neufeld (1996) shock 
model, indicates that at least 2 shock components are re- 
quired. One component at 15-20 km s _1 matches well the 
observed column densities in the v=0, J~3 to 7 levels, 
whilst a second component, at 35 km s _1 , is required to 
match higher J levels and the ro- vibrational lines (Wright 



1999). Harwit et al. (1998) used a shock velocity of 37 km 
s _1 to match their 8 LWS water line emission detections 
from 72 to 125 /um. These data were taken at a slightly 
different position, centred on BN rather than IRc2 and off- 
set 8.2" to the north. The SWS results may therefore be 
more sensitive to the low- velocity plateau gas identified by 
several authors (e.g. Gcnzcl & Stutzki 1989), whereas the 
larger LWS beam will include a significant contribution 
from the high-velocity gas. We also note that water line 
emission in the Kaufman & Neufeld (1996) shock models is 
for many lines insensitive to the shock velocity, especially 
when the velocity is > 15-20 km s _1 . This is presumably 
because in such cases the water abundance has reached 
its maximum value and that the gas temperature is well 
above the upper levels involved. Therefore, the emission 
depends more on density and not temperature (which in 
turn is highly sensitive to the shock velocity). Even so, 
future attempts to match both the Harwit et al. lines and 
our detections with a slower shock, but including radiative 
excitation, will be worthwhile. 



A long standing problem in shock research is the con- 
tribution of water to the total gas cooling behind the 
shock, and it is pertinent to mention it here. Unfortu- 
nately, the situation in Orion IRc2, i.e. a high water abun- 
dance (see Sect. 4.3) coupled with an intense mid- and 
far-infrared radiation field, has not previously been con- 
sidered in shock models. However, the detection of all 
our transitions in absorption raises the interesting pos- 
sibility that at least these transitions are not cooling the 
gas, but rather heating it. Such heating results from the 
water molecules absorbing the radiation and thereby be- 
ing raised to excited states. The excess energy may then 
be imparted to the gas through collisional de-excitation 
with other molecules, e.g. H 2 , as in the scenario proposed 
by Takahashi et al. (1983). This would of course be in 
competition with radiative de-excitation (i.e. cooling), and 
is therefore density dependent. Further, Takahashi et al. 
(1983) and Neufeld et al. (1995) state that a necessary 
condition for H 2 heating is that the dust temperature 
is greater than the gas temperature. A full scale shock 
model, taking into account all transitions, would be re- 
quired to determine if there was net heating or cooling. 
We merely mention it here as an interesting aside, and di- 
rect the reader to the paper by Harwit et al. (1998) for a 
full discussion of the water cooling contribution in Orion. 
Briefly though, from their larger beam LWS observations 
they find that H 2 and CO contribute similar amounts 
to the overall cooling, but both are only about a tenth of 
the total H 2 cooling. On the other hand, Saraceno et al. 
(1999), in a study of a sample of shock sources, find that 
CO cooling typically dominates over H 2 by a factor of 
several. 
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4-3. H 2 abundance in Orion 

In order to determine an H2O abundance directly from 
our data, information on the H 2 column density is needed. 
As noted above, from the H 2 lines detected elsewhere in 
our spectra (van Dishoeck et al. 1998), the shock has 
been characterized using the Kaufman & Neufeld (1996) 
C-shock models, in a manner similar to that applied by 
Wright et al. (1996) to Cepheus A. Further details can 
be found in Wright (1999); in short, the H 2 pure rota- 
tional lines originating from J=3 to J=7 indicate a col- 
umn A(H 2 ) = (1.2±0.2) x 10 21 cm" 2 of warm (T ex =700 K) 
gas, and the inferred shock velocity of 15-20 km s _1 is sim- 
ilar to the 18 km s _1 expansion velocity deduced from our 
FP observations, but lower than that obtained from mod- 
cling the H 2 vibration-rotation lines (35 km s" 1 ). How- 
ever, supporting evidence for the co-location of the H 2 
and shocked H 2 comes from Scoville et al. (1982), who find 
that at and around the position of IRc2 the peak emission 
of the 1-0 S(l) line is at vlsr=— 6 to —16 km s _1 , sim- 
ilar to the velocity of maximum H 2 absorption in our 
spectra. 

Using the above H 2 column density, and that of H 2 
calculated from the generalised curve-of-growth analysis, 
the inferred water abundance would be 1.25 x 10~ 3 . How- 
ever, the inferred column density from the high- J (i.e. 3- 
7) lines is a lower limit on the total column density, since 
the bulk of the H 2 is likely to be cooler than T ex =700 K. 
Such cooler gas would be probed by the 28.2 /zm H 2 0-0 
S(0) line, but which unfortunately has not been detected 
toward IRc2, due to the extremely high continuum emis- 
sion. However, a lower limit implies that T ex > 110 K, 
calculated from the populations in J=3 and 2, inferred 
from the observed 0-0 S(l) and upper limit S(0) lines re- 
spectively. On the other hand, the S(0) line is detected 
at about the 3<r level toward the shock "peak 2", ~ 30" 
to the south-east of IRc2, and in this case the excitation 
temperature between the J=3 and J=2 levels is ~ 150 K 
(Wright 1999). 

This temperature of ~150 K is likely to be indicative 
of the kinetic temperature, Tkin, at this position, since 
the density of >10 5 cm" 3 (e.g. Genzel & Stutzki 1989) 
is above the critical density of the H 2 0-0 S(0) and S(l) 
transitions (e.g. Le Bourlot et al. 1999). Although this 
estimate of Tkin is at a different position than where the 
water absorption line detections were made, Wright (1999) 
shows that in at least 3 shock sources, including Orion, the 
excitation temperature determined from pure rotational 
H 2 lines ( J up =3 7) is quite invariant across the sources, as 
is the column density. Other estimates of Tkin in the low- 
velocity plateau, tabulated by Genzel & Stutzki (1989), 
range from 100-500 K. 

If we assume that Tkin is between 100-500 K, then Ta- 
ble 2 presents the expected total H 2 column density based 
on the observed surface brightness of the 0-0 S(l) line to- 
ward IRc2 of 1.5xl0" 10 W cur 2 sr" 1 . We believe that 



Table 2. Inferred total H 2 column densities and water 
abundances 



Tahiti 


iV(H 2 ) a 

/ — 9 \ 

(cm 2 ) 


Ar(H 2 0)/iV(H 2 ) b 


(Kelvin) 




100 


1.1 xlO 24 


1.4xl0" 6 


150 


5.7xl0 22 


2.6xl0~ 5 


200 


1.4xl0 22 


l.lxlO" 4 


230 


8.1xl0 21 


1.9xl0" 4 


250 


6.2xl0 21 


2.4xl0" 4 


280 


4.4xl0 21 


3.4xl0" 4 


300 


3.7xl0 21 


4.1xl0" 4 


320 


3.2xl0 21 


4.7xl0" 4 


350 


2.6xl0 21 


5.8xl0" 4 


400 


2.1xl0 21 


7.1xl0" 4 


450 


1.7xl0 21 


8.8xl0" 4 


500 


1.5xl0 21 


l.OxlO" 3 



a Assuming that the observed H 2 0-0 S(l) 17.0348 /im line, 
with surface brightness 1.5xl0 -10 W cm" 2 sr _1 , arises 
from gas with kinetic temperature T kin . 

b Using a total water column density of 1.5x 10 18 cm" 2 , and 
assuming the H 2 and H 2 are co-located. 



100 K is too low to be applicable to the shocked plateau 
gas, since it is lower than Tkin we find at "peak 2" , and in 
any case would imply a column density an order of magni- 
tude larger than the value of 10 23 cm" 2 quoted by Genzel 
& Stutzki (1989). Also shown in Table 2 is the inferred wa- 
ter abundance, A(H 2 0)/A^(H 2 ), using the observed water 
column density of 1.5 x 10 18 cm" 2 , and assuming the water 
and molecular hydrogen to be co-located. Taking the max- 
imum possible water abundance to be 5 x 10" 4 (anything 
higher would violate the constraint imposed by the oxy- 
gen abundance itself of [O]/[H]=3.19xl0" 4 from Meyer 
et al. 1998), then the water abundance is in the range of 
2.6-50xl0" 5 . 

Let us instead assume a priori that our water detec- 
tions arise only in that portion of the shock where the 
temperature is high enough for the water creation reac- 
tions to rapidly proceed. Several authors (e.g. Hollcnbach 
& McKee 1979, Neufeld et al. 1995, Ceccarelli et al. 1997 
and Charnley 1997) have shown that water begins to be 
quickly formed when the temperature exceeds about 230 
K, and has reached its maximum abundance by about 
320 K. In that case our inferred water abundance is in the 
range 2 5xl0" 4 , in remarkable agreement with the result 
of Harwit et al. (1998) of 3.5-5 xlO" 4 considering the dif- 
ferent methods used. Further, it is interesting to note that 
Watson et al. (1985) find a H 2 column density of 3 x 10 21 
cm" 2 , derived from several observed far-infrared CO lines 
in a 44" beam, implying a H 2 abundance of 5xl0" 4 . 
Watson et al. however infer a kinetic temperature of 750 
K, close to the excitation temperature derived from our 
H 2 data (Wright 1999), whereas from Table 2 their col- 
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umn density would instead imply a temperature of only 
~330 K. 

Alternatively if we simply use a H 2 column density 
of 10 23 cm~ 2 for the plateau gas, as given by Genzel & 
Stutzki (1989), then the water abundance is 1.5 x 1CP 5 , 
in better agreement with that of other shocks (e.g. Liseau 
et al. 1996, Ceccarelli et al. 1998, Spinoglio et al. 1999), 
and suggesting that a substantial fraction of the gas-phase 
oxygen may still be in another form, most likely gas-phase 
atomic O I. However, this is probably a strict lower limit, 
since such a large column density probably includes gas 
at a temperature lower than 230 K, as well as gas behind 
the infrared continuum source. 

An H 2 abundance at the high end of the range quoted 
above, i.e. 2-5 xl0~ 4 , is favored if we assume that water is 
only efficiently created in regions where the temperature 
is greater than about 230 K. This water abundance is con- 
sistent with all of the gas-phase oxygen not locked up in 
CO being incorporated into H2O. However, H2O ice evap- 
orates at temperatures above about 90 K, so lower temper- 
atures, and hence a lower abundance, cannot be excluded. 
The major uncertainty on our abundance value is the total 
hydrogen column density. As mentioned above, our result 
is in good agreement with Harwit et al. (1998), but also 
in agreement with Cernicharo et al. (1999a, b) from LWS 
and ground based observations respectively. From radia- 
tive transfer modelling of their water detections, they find 
toward the Orion plateau a water abundance of order 1— 
2xl0 -4 , at most a factor of a few lower than what we find. 
Besides Orion, there is only one other source so far ob- 
served by ISO with such a high water abundance, namely 
L1448-mm by Nisini et al. (1999), with an abundance of 
order 5 x 10~ 4 . 

A similar treatment for the OH lines detected in our 
spectrum yields a column density of approximately 5 x 10 16 
cm~ 2 , after fixing T ex to be the same as that for water. 
The implied OH abundance is then (8.8 — 156) x 10 -7 for 
T fei „=150-320 K, or (6.2 - 15.6) x 10~ 6 for 230-320 K. 

5. Conclusions 

Through use of the ISO-SWS in its grating mode, 19 pure 
rotational absorption lines of water have been detected 
toward Orion IRc2 for the first time. Fabry-Perot spectra 
of 5 lines reveal that the water is located in an outflow 
expanding at a velocity of 18 km s _1 . The strong mid- 
infrared continuum toward IRc2 plays a dominant role 
in the excitation of the molecule and the line formation, 
which can be modeled using a simple, generalised curve- 
of-growth technique. This yields a total water column den- 
sity of order 1.5xl0 18 cm~ 2 and excitation temperature 
of 72 K, similar to the dust continuum colour tempera- 
ture. Both derived quantities have a maximum uncertainty 
of about 20%. The data provide support for large abun- 
dances of H2O in the outflows of massive stars. Simulta- 
neous analysis of the complete ISO-SWS and LWS data 



set on H2O, OH and CO may provide further informa- 
tion on the abundance and excitation of these molecules 
in the various physical components within the complex 
Orion environment. 
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